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ABSTRACT 

It is shown how to construct exactly gauge-invariant ^'-matrix elements for pro¬ 
cesses involving unstable gauge particles such as the W and bosons. The re¬ 
sults are applied to derive a physically meaningful expression for the cross-section 
cr(e+e“ ^ Z^Z^) and thereby provide a solution to the long-standing problem 
of the unstable particle. The problem of maintaining QED gauge-invariance in 
the process e+e“ ^ VeG~ud is examined. 

1. Introduction 

Pair production of W bosons has now been observed at LEP and both fU’s and 
Z°’s will be copiously produced at upcoming colliders. Yet both the W and are 
unstable with quite sizable widths that cannot be treated using the standard methods 
of S'-matrix theory because of the classic problem of the unstable particieS. Roughly 
stated the problem is as follows: S'-matrix theory deals with asymptotic in-states 
and out-states that propagate from and to infinity. Unstable particles cannot exist 
as asymptotic states because they decay a finite distance from the interaction region. 
Indeed it is knownil that the S'-matrix is unitary on the Hilbert space spanned by 
stable particle states and hence there is not even any room to accommodate unstable 
particles as external states. How can one use the S'-matrix to calculate, say, the 
production cross-section for an unstable particle when it cannot exist as an asymptotic 
state? 

In order to account for the finite width of the W and one normally performs 
a Dyson summation of the one particle irreducible (IPI) self-energy diagrams that 
effectively replaces the tree-level propagator by a dressed propagator. However the 
IPI self-energy is gauge-dependent starting at 0(a) and hence dressed propagator is 
gauge-dependent at all orders in a. If the dressed propagator is used in a hnite-order 
calculation the resulting matrix element will also be gauge-dependent at some order. 
This gauge-dependence should be viewed as an indicator that the approximation 
scheme being used is inconsistent or does not represent a physical observable and is 
a latter-day manifestation of the problem of the unstable particle. 

2. The Problem of the Unstable Particle 

In this section we obtain a physical and therefore gauge-invariant property that 


may be used to identify an unstable particle without requiring that it appear in the 
hnal state. 

Recall that the coordinate space dressed propagator for a scalar particle has an 
integral representation 


A(a;' 


d^k 

{271^ k'^ — TV? — n(/c2) + ie 


( 1 ) 


The integrand has a pole at k"^ = Sp where Sp is a solution of the equation s — m? + 
n(s) = 0 and an analytic function, F(s), via the relation s—m^+n(s) = {s—Sp)/F{s). 
The dressed propagator can then be written as 
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that separates resonant and non-resonant pieces. Performing the ko integration on 
the resonant part gives 
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where ko = \jk^ + Sp. The non-resonant term contributes only for t = t' and so 
represents a contact interaction. The resonant part spits into two terms that con¬ 
tribute when t > t' or t < t' and therefore connects points x and x' that are separated 
by a hnite distance in space-time. When a similar analysis is applied to a physical 
matrix element one concludes that the resonant part describes a process in which 
there is a hnite space-time separation between the initial-state vertex, Vi, and the 
hnal-state vertex, Vf. In other words, the resonant term describes the hnite propaga¬ 
tion of a physical unstable particle. The non-resonant background represents prompt 
production of the hnal state. As these two possibilities are, in principle, physically 
distinguishable, they must be separately gauge-invariant. 

We can thus use hnite propagation as a tag for identifying unstable particles 
without requiring that they appear in the hnal state. This is, after all, the way 6- 
quarks are identihed in vertex detectors. A production cross-section for an unstable 
particle is obtained by extracting the resonant part of the matrix element for a process 
containing that particle in an intermediate state and summing over all possible decay 
modes. 


3. The process e’''e — 

In this section we will calculate the cross-section for e’''e“ — Z^Z^. This is of both 
theoretical and practical importance. On the theoretical side it represents an example 



of a calculation of the production cross-section for unstable particles. On the practical 
side, at high energies will be a dominant source of fermion pairs (/i/i) 

and (/2/2) due to its double resonant enhancement and hence a{e'^e~ is 

an excellent approximation to the cross-section for 4 -fermion pair production. If the 
experimental situation warrants it, background terms can also be included without 
difficulty. A more detailed account can be found in ref. ||. 

In the case of //, dealt with in ref.s the invariant mass squared of 

the s, is hxed by the momenta of the incoming For the process e’''e“ —>■ 

Z^Z^ the invariant mass of the produced Z°’s is not constant and must be somehow 
included in phase space integrations. It is not immediately clear how to do this and 
without further guidance from S'-matrix theory there would seem to be considerable 
flexibility in how to proceed. A new ingredient is required and that is to realize that 
an expression for an S'-matrix element can always be divided into a part that is a 
Lorentz-invariant function of the kinematic invariants of the problem and Lorentz- 
covariant objects, such as etc. The latter are known as standard covananfsiffli. 
It is the Lorentz invariant part that satishes the requirements of analytic S-matrix 
theory and from which the resonant and non-resonant background parts are extracted 
while the Lorentz covariant part is untouched. 

To calculate a{e^e~ —>• Z^Z^), we begin by constructing the cross-section for 
Z^Z^ (/i/i)(/2/2), and will eventually sum over all fermion species. The 


e' e 


squared invariant masses of the /i/i and 72/2 pairs will be denoted pf and p^- 

To obtain the piece of the matrix element that corresponds to finite propagation 
of both Z°’s we extract the leading term in a Laurent expansion in pj and of th® 
analytic Lorentz-invariant part of matrix element leaving the Lorentz-covariant part 
untouched. This is the doubly-resonant term and is given by 


i 

^ ['»/i7^(f4/Jgp)7L + Vzfj,{sp)-fR)vf^] (3) 

Pi 

^ [ff/27'^(f4/Jgp)7L + Vzf^{sp)'yR)vf^] 

P2 ^p 

where are Lorentz covariant tensors that span the tensor structure of the matrix 
element and 7^, are the usual helicity projection operators. The Mj, fl^^ and Vzf 
are Lorentz scalars that are analytic functions of the independent kinematic Lorentz 
invariants of the problem. Fzz is defined in the way described following eq.( 2 ). It 
should be emphasized that eq.(^) is the exact form of the doubly-resonant matrix 
element to all orders in perturbation theory that we will now specialize to leading 
order. 

In the procedure described above, one starts by extracting the resonant term 
in a scattering amplitude by Laurent expansion about the exact pole position Sp 




and then specializes to lower orders by farther expanding abont the renormalized 
mass. Other anthors0’0 have attempted to apply the techniqnes described above 
by first expanding abont the renormalized mass and then adding a finite width in 
the denominator of the resonant part by hand. That procednre cannot be justified 
and leads to problems when one treats processes like e+e“ —>• It gives rise 

to spurious threshold singularities or complex scattering angles due the production 
threshold’s branch point being incorrectly located on the real axis. 

It is important to note that Laurent expansion refers to how the resonant term 
is identihed. The aim is only to separate M. into gauge-invariant resonant and back¬ 
ground pieces and there is no necessity to perform a Laurent expansion beyond its 
leading term although this can provide a useful way for parameterizing electroweak 
datall3. Although in what follows the background terms are systematically dropped, 
they can be kept if the experimental situation should require it. 

In lowest order eq.(|D becomes, up to overall multiplicative factors. 
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where Mi = t = 7 j.(|ie- -^ 2 ) 11 ,, M 2 = u ^ and the hnal 

state vertex corrections take the form Vzf^ = iejS^jL and Vzf^ = The left- 

and right-handed couplings of the to a fermion / are 
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Squaring the matrix element and integrating over the hnal state momenta for 
hxed Pi and gives 
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where the sum is over fermion species. Note that —>■ 6{p‘^—M^)0{po) as Im(sp) —> 

0 which is the result obtained by cutting a free propagator. The variables s, t, u, 











pI and P 2 in eq.(||) arise from prodncts of standard covariants and external wave 
fnnctions and therefore take real valnes dictated by the kinematics. 

Integrating over t, p\ and p^ leads to 


where 
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4. The Process e’''e — h>e& ud 

Massless particles, snch as the photon, are always problematic for analytic S- 
matrix analyses. In this section the process e’''e“ —> iJf,e~ud will be considered with 
particnlar emphasis on how take into acconnt the hnite width of the W boson while 
at the same time maintaining QED gange invariance. The dominant contribntion 
to this process is expected to come from physical W~^ prodnction in which e’''e“ —> 
Uee~W~^ —3- Uee~ud via diagrams of the type shown in Fig.l. Assnme for the moment 
that the W is stable and therefore has the nsnal propagator (p2 + M^)-F It is 
expected that the matrix element for this process shonld satisfy the nsnal condition 
for QED gange-invariance 

= 0 (7) 

where is the matrix element with the external electromagnetic cnrrent and photon 
propagator ampntated and is the momentnm in the photon propagator. When this 
test is applied to the diagrams of Fig.l. one hnds that (|^ does not hold. The reason 
is that there exist other Feynman diagrams that prodnce the same hnal state and 
only when these are inclnded is the condition (|^ correctly satis£ed0. However the 
additional diagrams that mnst be inclnded do not contain internal IT’s that conld go 
on-mass-shell and are therefore expected to be snbdominant. 

It might be thonght that gange dependence in snbdominant terms of no concern 
bnt for the process in qnestion it plays a crncial role. It can be shown0 that as 0 
the resnlting cross-section behaves as ~ 1/g^ provided (0) is satished and ~ l/g*^ when 
it is not. Withont the maintenance of QED gange invariance the cross-section then 
blows np too strongly as g^ —> 0. 

The incinsion of snbdominant diagrams to cancel the gange-dependence in the 
dominant contribntion contrasts sharply with what was done in section 3. The Lan- 
rent expansion extracts an exactly gange-invariant resonant term and eliminates from 








Figure 1: Tree-level diagrams contributing to the resonant part of the process e’''e 
Ve&~ud. (courtesy of the authors of ref.|l^ 


it any non-resonant (subdominant) pieces. These pieces that are removed may be 
ultimately combined with other contributions to yield an overall gauge-invariant non¬ 
resonant contribution if required. 

Now consider the case when the hnite width of the W boson in the diagrams of 
Fig.l is to be taken into account. If the procedure described in sections 2 and 3 is 
applied to the process —>■ —>■ Uee~ud to extract the dominant resonant 

piece then the resulting matrix element does not satisfy the condition (J^) and the 
resulting cross-section blows up too strongly as —*• 0. A number of solutions 
have been proposed and their relative merits examined in ref. 


14. One that has 


received much attention is that of Baur and Zeppenfeldlil. Their solution is to add 
the imaginary parts of certain higher-order diagrams, specihcally the 1-loop fermionic 
correction to the W^W~'y vertex, which they show exactly removes the QED gauge 
dependence. This method has the practical disadvantage that it requires an additional 
calculation that is signihcantly more complicated to obtain than the rest of the tree- 
level amplitude. Moreover once it is calculated it is immediately largely canceled by 
the gauge-dependent terms in the resonant part. 

In this section we describe an alternative more economical approach that lies 
more in line with the strategies adopted in section 3. Instead of seeking additional 
higher-order gauge-dependent corrections to cancel the QED gauge dependence, we 
eliminate it by means of a projection operator under which the exact matrix element 
is known to be invariant. 

Suppose the approach of section 3 has been used to obtain the leading resonant 
contribution to the process —>■ —>• Ve^~ud. After amputating the ex¬ 

ternal electromagnetic current and photon propagator a Green’s function, is ob¬ 
tained that can be written in terms of a set of standard covariants, /f, and associated 
Lorentz scalar functions, Ai{p\,p 2 ^ 




( 8 ) 


where the arguments pi, p 2 , etc. are 4-momenta that appear in the problem. The /f, 
of course, must transform in the same way that does. 



The gauge condition on Green’s functions qfj,G^ = 0 means that the Ai are not 
linearly independent since they satisfy the condition ' h)Ai = 0. 

Gauge-invariant Green’s functions can be obtained using the methods of Bardeen 
and Tun^. Suppose that we have some exact Green’s function expressed in terms 
of standard covariants and Lorentz scalars as in eq.(|). Since = 0 the Green’s 
function is invariant under of the action of the operator 
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where p is some conveniently chosen momentum. The spinor-Lorentz tensor structure 
of the Green’s function is spanned by the set That is 


g'‘(pi,P2, ■■■) = Y, (^^*!'(pi.P2. ■■■)) Adi.pi 
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For an approximate Green’s function that does not satisfy the QED gauge-condition 
the operator serves as a projection operator onto From the it 

is always possible to construct a new basis that is free of kinematic singularities that 
might have occurred when, for example, {p-q) = 0. In practice, however, it is simpler 
to just choose the vector p such that {p ■ q) is not small in the kinematic region of 
interest. 

In general consider a Green’s function calculated using some incomplete ex¬ 
pansion up to some order. It will consist of two parts 


+ Gj' (11) 

where Gg is a consistent gauge-invariant contribution correct to given order of the 
calculation. Gj* is a spurious higher-order gauge-dependent correction. Because Gg 
satisfies the gauge condition, g^^Gg = 0, it is invariant under the action of the operator 
/)(. Thus at a given order we may make the replacement 

Gm ^ jMG^ = G^ + J^Gj'. 

Only Gi is affected by the projection operator but since it is of higher order this 
is of no concern. The important point is that it now satisfies the gauge condition 
g^(/^Gj') = 0 and is prevented from causing too strong a numerical blow up as 
g^ ^ 0. 

Now, of course, the result depends on the choice of the vector, p, which when 
judiciously done can be used to significantly simplify calculations. The dependence 
on p, as it appears only in higher-orders is analogous to renormalization scheme 
dependence, or the p-dependence in the MS scheme, that is an inescapable feature 
of any perturbative expansion. As with scheme dependence, we can be sure that it 
will become ever weaker as higher orders included since the exact Green’s function is 
invariant under the action of 



Details of applying the projection operator to the process > 

Vee~ud may be found in ref. but may be summarized as follows: 

The methods of section 3 are used to obtain the resonant part of the matrix element 
M. which takes the form 

M = JjPiG^ (12) 

where Jj is the electromagnetic current of the electron and is the photon propa¬ 
gator. It is simply a matter of inserting the projection operator, (^, to obtain 

M = JlP^ir^G^). (13) 

The terms in parentheses satisfy the gauge-condition (|^ and therefore the cross- 
section formed from M. behaves ~ 1/g^ as —>■ 0. Terms eliminated by the projection 
operator must cancel in the exact matrix element. 
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